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ABSTRACT 

Aims. We build a simple analytical model for the bias of dark matter halos that applies to objects defined by an arbitrary 
density threshold, 200 < 5* < 1600, and that provides accurate predictions from low-mass to high-mass halos. 
Methods. We point out that it is possible to build simple and efficient models, with no free parameter for the halo 
bias, by using integral constraints that govern the behavior of low-mass and typical halos, whereas the properties of 
rare massive halos are derived through explicit asymptotic approaches. We also describe how to take into account the 
impact of halo motions on their bias, using their linear displacement field. 

Results. We obtain a good agreement with numerical simulations for the halo mass functions and large-scale bias at 
redshifts < z < 2.5, for halos defined by a nonlinear density threshold 200 < S, < 1600. We also evaluate the impact 
on the halo bias of two common approximations, i) neglecting halo motions, and ii) linearizing the halo two-point 
correlation. 

Key words. Cosmology: large-scale structure of Universe; gravitation; Methods: analytical 



1. Introduction 

A fundamental test of cosmological models is provided 
by the distribution of nonlinear virialized objects, such as 
galaxies or clusters of galaxies. First, this allows to check 
that nonlinear objects form through the amplification by 
gravitational instability o f small (nearly Gaussian) primor- 
dial density fluctuations (|Peebles!ll980[ ). Second, quantita- 
tive comparisons between theoretical predictions and ob- 
servations allow to derive constrain ts on cosmological pa- 
rameters (jZheng fe Weinberg! I2007T) . Indeed, observations 
show that galaxies and clusters do not follow a Poisson dis- 
tribution but show signi fi cant large-sca l e corre lations (e.g., 
iMcCracken et all (|2008l ): iPadilla et all (|2004f l). In partic- 



ular, at large scales their two-point correlation function 
is roughly proportional to the underlying matter correla- 
tion, with a multiplicative factor 6 2 , called the bias, that is 
greater for more massive objects. 

Following the spirit of the Press-Schechter picture 
(|Press fc Schechterl I1974D , where nonlinear virialized ob- 
jects are identified with h igh-density flu ctuations in the ini- 
tial (linear) density field, iKaiserl (|1984[ ) showed how such a 
behavior for the halo correlation naturally arises. Similar re- 
sults are obtained for the clustering of densit y peaks instead 
of ov erdensities above a given threshold (jBardeen et al.l 
|1986() . Indeed, distant high-density regions are correlated 
through their common longwavelength modes of the lin- 
ear density field, and the effect is greater for more massive 
and extreme objects. A somewhat simpler derivation can 
be obtained through the peak-backgro und split argument 
(|Cole fc Kaiserlll989t [Bond et al.lll99lH . This rests on the 
same physics, which is that within a large region character- 
ized by a positive mean density contrast local density peaks 
beyond a given density threshold are more frequent than in 
the mean, which yields a positive correlation between rare 
massive halos and the larger-scale matter density field. 



Along these lines, a simple analytical model that com- 
pares reas onably well w i th num erical simulations was pre- 
sented by IMo fc White] (119961) . wh ile a better agreement 
was obtained by Sheth fc Tormenl (|1999|) . at the cost of 
adding a few parameters fitted to the simulations (through 
the halo mass function). More complex models, that also 
apply to redshift space and take into account sublead- 
ing scale-d e pende nt terms can be found for in stance in 
iDesiacauesl (|2008f k iDesiacaues fc Shethl (|2010h . On the 
other hand, various fitting formulas have been proposed 



from measures in numerical simulations J Hamana et al 
| 200lblPiUepich. et al l20TfjHTinker et al l201fj|Manera et al 



2010). In particular. lTinker et al.l (|2010l ) have recently stud- 
ied the dependence of the halo bias on the density contrast 
(5* used to define the halos. In the paper we present a sim- 
ple analytic al model for th e halo bias that extends a pre- 
vious work (Valageas 2009) to arbitrary density thresholds 
in the range 200 < <5» < 1600. This is of practical interest 
as the definition of halos can vary among authors, and it 
can happen that observations only extend to smaller radii 
than the usual virial radius, which corresponds to objects 
defined by higher density t hresholds. In addition, we also 
improve on IValageasI (2009) by simplifying somewhat the 
model, especially the treatment of halo motions, while ob- 
taining a better accuracy for low-mass halos. This should 
allow an easier generalization to more complex cases, such 
as non-Gaussian initial conditions. This simple model also 
allows us to evaluate the inaccuracies implied by two ap- 
proximations that are frequently used, that is, neglecting 
halo motions and linearizing the halo correlation function 
over the matter correlation function. 

In Sect. [5] we describe our model for the bias of halos de- 
fined by a nonlinear densi ty contrast 5* = 200. We explain 
how the standard model oflKaiserl (19841) for the correlation 
of rare massive halos is modified once we take into account 
the motion of halos. Then, we point out that normaliza- 
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tion conditions can be very useful to constrain the bias of 
intermediate and low mass halos, which allows to build a 
simple and efficient model. Next, we show that we obtain a 
good agreement with numerical simulations. In Sect. |3]we 
extend this model for the bias, as well as for the halo mass 
function, to the case of halos defined by a larger nonlinear 
density contrast, and we compare our results with numer- 
ical simulations for 200 < S* < 1600. Finally, we estimate 
the importance of halo motions and of the nonlinearity of 
the bias in Sects. Hand 02 and we conclude in Sect. [5] 



2. Bias of dark matter halos 



Following iKaiserl (jl984h we estimate the two-point correla- 
tion function of dark matter halos, whence their bias, from 
the probability to reach a linear density threshold 5l* (as- 
sociated with the formati on of viria l ized h alos of nonlinear 
density contrast S*). As in Valageas (2009]), we improve this 
model by taking into account the motion of halos, associ- 
ated with the change from Lagrangian to Eulerian space 
(i.e., from the linear density field to the actual nonlinear 
de nsity field). How ever, we simplify the prescription used 
m IValageasI ipOOl as we no longer use a spherical collapse 
model to estimate these displacements but use the initial 
momenta of both halos (i.e., the linear displacement field). 
In addition, we add a further ingredient to the study of 
Valageas (2009), as we explicitly enforce the normalization 
to unity of the halo bias (when integrated over all halo 
masses). We describe below this simple model. 

2.1. Lagrangian space 

We first consider the two-point correlation function of halos 
in Lagrangian space, that is within the linear density field 
<5i(q), where we identify future halos of Eulerian radius r 
and nonlinear density contrast <5* (with typically S* ~ 200) 
with spherical regions of Lagrangian radius q and linear 
density contrast Sl*. Thanks to the conservation of mass, 
the Lagrangian and Eulerian properties of halos of mass M 
are related by 



M 



-in 



pq 3 , 



(l + 5*)r 3 and <5* = F[6 L +), (1) 



where p i s the mean matter density. As pointed out in 
IValageasI (|2009h . the function T describes the spherical 
collapse dynamics, so that in a ACDM cosmology with 
fi ro = 0.27, 0,^. — 0.73, we have 5l* — 1-59 at z = for 
(5* = 200, instead of the usual value S c ~ 1.675 associated 
with full collapse to a point, that is with 6* — oo. 

Defining halos by the linear threshold Sl* = J-~ 1 (6*), 
rather than by the full collapse value S c , is closer to obser- 
vational and numerical procedures, since (in the best cases) 
one defines "halos" in galaxy or cluster surveys, and in nu- 
merical simulations, by the radius and mass of overdensities 
within a given nonlinear density threshold d*. Moreover, 
this gives the freed om to chose d ifferent thresholds, such 
as ft = 200 o r 100 (|Valageasll2d09h . As discussed in Sect.3 
in IValageasI (|2009h . for massive (M ~ 10 15 /i 1 Mq) and 
rare halos, the choice (5* = 200 also roughly corresponds 
to the separation between outer shells, dominated by ra- 
dial accretion, and inner shells, with a significant trans- 
verse velocity dispersion, that have experienced shell cross- 
ing. The case of typical halos (i.e. below the knee of the 



halo mass function) is more intricate as they do not show 
such a clear separation. This can also be seen in numerical 
simulations, such as Fig. 3 in lCuesta et al.l ([2008). At high 
redshift, where nonlinear objects have a smaller mass, this 
"virialization" radius shifts to higher density contrasts for 
ACDM cosmologies, because of the change of slope of the 
linear matter power spectrum with s cale (e.g., ft ~ 5 00 for 
M ~ lO 11 /i- 1 M , as seen in Fig.5 in IValageasI (|2009f) ). 

Then, defining the halo correlation as the fractional ex- 
cess of halo pairs (|Kaiserlll984l: lPeeblesHl98(l . we write in 
Lagrangian space 

n L (Mi, M 2 ; si, s 2 )dA/idM 2 dsids2 = 

n L {M 1 )n L {M 2 )[l + £ L {M X ,M 2 ; s)]dM 1 dM 2 ds 1 ds 2 . (2) 

Here and in the following we use the letter s for the position 
of halos in Lagrangian space to avoid confusion with their 
Lagrangian radius q, and we introduced the Lagrangian dis- 
tance, s = |s 2 — si I , between both objects. We also note 
with a subscript L the quantities associated with the lin- 
ear fie lds o r the Lagrangian space. Then, following IKaiserl 
(|l984f) and IValagcas (2009), we obtain the fractional ex- 
cess of halo pairs from the bivariate density distribution 
T- > l{^lii^L2) over the two spheres of radii q\ and q 2 , 

Vl(Sli,S L 2) 



1 + £l {Mi , M 2 ; s) 



<7i<7 2 



exp 



Pl(Sli)Pl(S L2 ) 



'12) 



(3) 



(4) 



where we considered a single population, that is halos de- 
fined by the same density thresholds 5l* and (5* = J~(6l*)- 
Here we assumed Gaussian initial conditions and we intro- 
duced the cross-correlation of the smoothed linear density 
contrast at scales q\ and q 2 , at positions si and s 2 , 



,0) = (<5 Lgi (si)(5 L?2 (si +s)} 



= 4tt / dkk 2 P L {k)W{k qi )W(kq 2 ) 
Jo 



sin(fcs) 
ks 



(5) 



where W{kq) is the Fourier transform of the top-hat win- 
dow of radius q, 



W{kq) 



dq ; 



kq 



t sm(kq) — kqcos(kq) 



(6) 



IV V {kqf 
and Ph{k) is the linear matter power spectrum, defined by 

S L (s) = J dke ik s 5 L (k), (7) 

(5 L (k 1 )5 i (k 2 )) = a D (k! + k 2 ) P L (h). (8) 

In particular, <j q = cr q>q (0) is the usual rms linear density 
contrast at scale q. In Eq.Q we also used the short-hand 
notation o\ = tf uqi (0), v 2 = ^ 2<q2 (0), and u\ 2 = (T^ >Qa (s), 
for the covariances over the two spherical Lagrangian re- 
gions. 

Of course, the prescription (HJ) only applies to the limit 
of rare events (i.e., large mass M) and large separations, 
where halos are isolated and almost spherical, so that we 
can neglect tidal effects as well as violent processes such 
as mergings. We shall come back to this point in Sect. 12.31 
below. The Eq.Q expresses the fact that if we have a rare 
overdensity at position si the probability to have a sec- 
ond overdensity at a nearby position s 2 is amplified, with 
respect to random locations, because of common longwave- 
length modes. 
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2.2. Eulerian space 

We must now convert the expression (j4]) into a model for the 
halo correlation in Eulerian space, that is within the actual 
nonlinear density field. Indeed, even halos that have not 
been destroyed or strongly modified by merging events have 
had time to move since the early times associated with the 
linear density field (which formally corresponds to z — > oo). 
In particular, the particles they are made of have moved 
even before the formation of the final halo. This leads to two 
effects with respect to the halo correlation: i) the Eulcrian- 
space distance x between the two halos is not identical to 
the Lagrangian-space distance s, ii) the halo number densi- 
ties are modified because the Lagrangian-Eulerian mapping 
modifies lo cal volumes (i.e ., dx ^ ds). These effects were es- 
timated in lValageai (|2009h from a spherical collapse model. 
Here we present a simpler method, which has the advan- 
tage of being more flexible and more closely related to the 
estimate (HJ), based on the linear density field. Moreover, it 
avoids the need to use a spherical dynamics approximation. 
As a further simplification, we also neglect the second ef- 
fect ii), due to the change of volume between Eulerian and 
Lagrangian spaces, because it is a subdominant prefactor. 

Thus, we simply estimate the displacement \P(M, s) of 
a halo of mass M and Lagra ngian positi o n s fro m the linear 
displacement field (see also iDesiacauesl (|2008h for a some- 
what different implementation in the case of density peaks). 
This gives for the Eulerian position of its center of mass, 



where we introduced the covariance 



x(M,s) = s + * L (M,s), 
with 

9 L (M,a)=i f dke iW s ^~5 L (k)W(kq). 



(9) 



(10) 



This coincides with the use of the Zeldovich approxima- 
tion (Zeldovich 1 970[) for the dynamics of the particles that 
make up the halo of mass M. However, since the halo scale 
q introduces a natural cutoff at high k (through the factor 
W(kq)), which for large masses and rare events is within 
the linear regime, the approximation should fare much 
better than what could be expected from an analysis of the 
Zeldovich dynamics for individual particle trajectories. For 
instance, the "truncated Zeldovich approximation" , where 
one suppresses the initial power at high k, provides a good 
description of large-scale clustering and yields a significant 
impr ovement over the unsmo othed Zeldo vich approxima- 
tion (jColes et alJll993tlMelott et al.lll994h . Physically, this 
means that small-scale nonlinear processes such as virial- 
ization within bound objects conserve momentum and do 
not affect the large-scale properties of the system. 

Then, for Gaussian initial conditions we can obtain 
the trivariate probability distribution Vl(8li, &L2, ^LV2\\)> 
where VE , 12 || is the linear relative displacement of the halos 
along their Lagrangian separation vector, 



M2|| 



[^L(M 2 ,S 2 ) - ¥ L (Mi,Bi)] • (S 2 - Si 
|s 2 - Si I 



(11) 



In particular, the mean conditional relative displacement, 
at fixed linear density contrasts 5li and 8l2 within the two 
spheres of radii q± and q 2 , reads as 



(*£12||)*z,i,«i 



S L1 (a^aj 2 ) + 5 L2 (aj-(7f 2 ) 



^11 



2^-2 



-(<5.L1*L12||) = -(£e,2*L12||), 



which writes as 



2 



with 



3 a qi,i2,s 



(13) 



(14) 



< 92S = 47r/ dkk 2 P L (k)W(k qi )W(kq 2 )W(ks), (15) 
Jo 

where again s = |s 2 — si|. We introduced a minus sign in 
the definition (| X3[) to avoid a minus sign in Eqs. (fl4l - p3|) . 
Although a qi g s is not guaranteed to be always positive, it 
is usually positive, which expresses through (| L3[) that in the 
mean two rare overdensities tend to move closer, as could 
be expected through their mutual gravitational attraction. 
By symmetry, the mean conditional relative displacement 
along the orthogonal directions vanishes, Li2±)s L1 ,8 L2 = 
0. Then, at lowest order the mean Eulerian-space distance 
x between both halos reads as 



(x(M 1 ,M 2 ;s)) = s - ^d L ,a 2 quq2iS - 



(16) 



l u 2 



'12 



where we considered a single population defined by the 
same density threshold Sl*- Neglecting the width of the 
distribution Vl(^li2\\)i whence of Vl{x), and in the limit 
of large separation, s — » oo, where & qiqoS — > 0, we can 
invert Eq. (ITB|) as 



s(Mx,M 2 ;x) ~ x 



1 + 



S 



'l u 2 



'12 



, (17) 



where a\ 2 is also evaluated at distance x. 

It is interesting to note that the express ion (fTT)) is very 
close to the one obtained in ValageaJ (|2009f) from a spheri- 
cal collapse model. The main difference is that this new re- 
sult involves the quantity a qi q2 s defined in Eq. (ll5[) . which 

contains the three windows W(kqi), W(kq 2 ), and W(ks), 
whereas the former result only involved quantities such as 
<Tg S and a q (s) that only contain two windows. The rea- 
son is that the model used in IValaeeasl (|2009l ) treated each 
halo as a test particle within the spherical gravitational 
field built by the other halo, whereas in the derivation of 
Eq. (|17p above we did not need to make this approximation 
and we simultaneously take into account the effect of both 
halos on the initial gravitational field. Nevertheless, it is 
reassuring that both models give similar results, as could 
be expected from qualitative arguments (for rare and well- 
separated halos). 

Equation (|17l) gives the "typical" Lagrangian-space dis- 
tance s between two halos that are observed in the non- 
linear density field at the Eulerian distance x. This takes 
care of the first effect i) associated with halo motions. As 
noticed above, the second effect ii) of halo displacements 
is to change volumes so that dx ^ ds. Then, by conser- 
vation of the number of pairs, the Lagrangian-space and 
Eulerian-space correlations at distances s and x are related 
by 



err a 



,(12) 



1"2 



'12 



[1 + e(Mi, M 2 ; ar)]dx - [1 + fr(Afi, M 2 ; s)]ds, 



(18) 
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and we could use again Eqs.© and ([T7)l to estimate the fac- 
tor |ds/dx|. However, for rare massive halos, of Lagrangian 
radius q, this is a subdominant effect, since from Eq. (|17[) 
we have |As/Aa;| ~ SL*o- 2 qs /a 2 , whereas the argument of 
the exponential Q behaves as 8 2 L ^a 2 2 l u 1 ^. Thus, the latter 
is larger than the former by a factor 8l*/o~ 2 , which is much 
greater than unity (and this is further amplified for very 
rare halos by the exponential). Therefore, to avoid intro- 
ducing unnecessary approximations and to make the model 
as simple as possible, we neglect this volume effect and we 
only keep the exponential behavior of the two-point corre- 
lation of rare halos. This gives for equal-mass halos, with 
Ml = M 2 = M and qi = q 2 = q, 

l + £(M,x) ~ e ^* CT L«/W+<,«)], (19) 
whereas the Lagrangian distance reads from Eq. (|17p as 



s(M,x) = x 1 + 



2 8l* o\^ x 

3(^+^,(3)) 



(20) 



Then, since at large distance the dark matter two-point 
correlation function is within the linear regime and reads 
as £(x) ~ <7q q (x), we obtain for the bias of halos of mass 
M, 



b 2 ,JM,x) 



!i(M,x) 



>)/K 2 (^+<>))] 



(21) 



(22) 



Here the subscript "r.e." stands for the "rare-event" limit 
to recall that the results obtained so far have been derived 
for very massive and rare halos, and do not apply to small 
objects. 

2.3. Normalization 

It is very difficult to derive an approximate model for the 
bias, and even the mass function, of small halos. Indeed, 
such objects associated with typical density fluctuations 
have been strongly distorted by tidal effects and have often 
undergone various merging events. Then, it is not possi- 
ble to identify in a simple manner the precursors of these 
halos in the initial linear density field. In this article we pro- 
pose the following simple prescription to bypass this prob- 
lem: we add to the asymptotic result (1221) a constant term 
with respect to the halo mass, bo(x), which we obtain from 
the normalization of the bias. More precisely, making the 
approximation that the bias factorizes (which however is 
not exact, even in the lar ge-mass regime, see Eq. (|2"2")l and 
Fig. 12 in lValageasI (|2009h ). as 



b 2 (Mi ,M 2 ;x)~b(M 1 ,x)b(M 2 ,x), 
we model the bias b(M, x) as 
b(M,x) = b T . e .{M,x) +b (x). 



(23) 



Here 6 r . e .(M, x) is given by Eg. (|2"2l (taking the square-root 
of this positive quantity), while bo(x) is set by the normal- 
ization (at fixed distance x) 



i j 
b(M,x)f(u) — = l 



whence 
b {x) = 1 - 



b r . e .(M,x)f(p) — 
o v 



where we used 



f{ v ) — 
v 



1. 



(26) 



(27) 



Here we introduced the scaled mass function of dark matter 
halos, that is, we write the halo mass function as 



n (M)dM=^/( i /)^ 



with v = 



>L* 



a(M) ' 



(28) 



where Sl* is the linear density threshold that defines these 
halos, given by Eq.((T]), and a(M) = a q . The well-know 
constraints ((27)) and (|25|) follow from the requirement that 
all the mass is accounted for by the mass function of halos 
(whatever the choice of Sl*), which leads to Eq. (|2"T]) . so that 
the density field on large scales (beyond the size of these 
halos) can be written in te rms of the halo mass function, 
as in the usual halo model (jCoorav fc Shethl[2002t) (the so- 
called "two-halo term"). Then, requiring that one recovers 
the dark matter two-point correlation function leads to the 
constraint (|2"5|) . provided the bias can be factorized as (|2"3"f 
(otherwise the constraint involves a bidimensional integral 
over & 2 (Ml,M 2 )). 

We do not claim here that the halo bias has a nonzero 
asymptote at low mass. In fact, the asymptotic behavior 
of the low-mass tail of the bias (and of the mass function 
itself) is largely unknown, but numerical simulations show 
that the dependence on mass flattens below v ~ 1 and 
seems roughly constant down to v ~ 0.3. Then, the pre- 
scription (|24p is intended to describe both this behavior 
(since from expression (|2"21 we can see that 6 r . e .(M) — > 
for M — > 0, or more precisely when a q — >• oo) and the 
constraint (1231) . Indeed, if the asymptotic behavior (f2"2"|) 
provides a good description for v > 1 we can expect the 
constraint (1231) to provide a good estimate for the (almost 
constant) value of b(M) in the regime v < 1. Moreover, 
it is always useful to make sure normalization constraints 
such as (f2"5")) are satisfied by the models. Indeed, this en- 
sures that the models are self-consistent and that absurd 
results will not be produced by the violation of basic inter- 
nal constraints. 

Another advantage of this simple prescription is that 
our model (|24[) of the halo bias has no specific free param- 
eter, apart from those already contained in the halo mass 
function (especially its low- mass tail). This is why we prefer 
to keep a simple constant term bo, instead of introducing for 
instance higher-order polynomials (over v or M) that would 
require some fitting over numerical simulations. This pro- 
vides a greater flexibility to the model, which can be used 
for a variety of cosmologies. 

Finally, for n umerical compu tations we use the mass 
function given in iValaeead (pool . 



( 24 ) f(u) = 0.502 [(0.6 v) 2 -" + (0.62 ^f 5 ] e^ 2 , 



(29) 



(25) 



which has been shown to agree with numerical simulations 
(for halos defined by (5* — 200). The exponential falloff, 

e~ v / 2 , where v is defined by Eos.([28| and ([T]), is consistent 
with the exponential term of Eq.Q. Indeed, both the 1- 
point distribution (i.e. the mass function) and the 2-point 
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Fig. 1. The halo bias b(M,x) as a function of cr(M), at distance sc = 50ft.~ 1 Mpc, for halos defined by the nonlinear 
density contrast <5» = 200. We plot our results at redshifts z = (left panel), z — 1.25 (middle panel), and z = 2.5 
(right panel). The black solid line, "<5l*"> is the theoretical prediction (j2"4"l) . the red dashed line, "s = a;", corresponds 
to the approximation where halo motions are neglected, and the gree n dot-dashed line, "S c " , corresponds to the use of 
5 C = 1.686 instead of Sl*- The blue dashed l i ne, "S MT" , is the fit from lSheth et al.l ((2001) , and the points are the results 
from numerical simulations of iTinker et al.l (|2010f ) . 



distribution (i.e. the halo correlation or bias) are obtained 
in the large-mass limit from spherical overdensities 5l* in 
the linear density field, with Sl* = J r_1 (<5*) and J* = 200. 
Note also that the normalization of the mass function (|2"9")l 
is not a free parameter since it is set by the constraint (|27p . 
This is why we prefer to use the mass function (j2"9")l , so that 
the model is fully self-consistent and large-mass tails do not 
involve free parameters. 

In particular, following ICole fc Kaiser! (|l989t) . apply- 
ing the peak-background split argument to Eq. (j2"9")l gives 
b ~ <5_l*/(t^, in the rare-event and large-distance limits. 
This agrees with the asymptotic behavior of Eq. (j2"2")) . ex- 
cept that Ea. (|2"2"]) also yields the prefactor a qvq (s)/<7ofl(x), 
which is different from unity. This expresses the facts that 
i) contributions to the correlation of objects of size q are 
damped for high wavenumbers, k 3> l/<?, and ii) halo mo- 
tions are nonzero and correlated, s ^ x. Of course, these 
two effects are neglected by the peak-background split ar- 
gument. 

2.4. Comparison with numerical simulations 

We now compare the model defined by Eqs. (f2"U|) . (|22[). 
(j2"4"]l . and (l2l)l) . with results from numerical simulations. 
We first consider in Fig. [T] the dependence on halo mass 
M of the large-scale bias, for halos defined by the non- 
linear density contrast <5* = 200 (whence 8l* — 1.59) at 
distance x — 50/i _1 Mpc. We show our results at redshifts 
z = 0,1.25, and 2.5. We can see that we o btain a good 
agreem ent with the numerical simulations of ITinker et al.l 
(12010ft . This was expected at high masses, where the ar- 
guments based on the clustering of rare overdensities in 
the lin e ar Ga ussian density field apply (as introduced by 
iKaiserl (|l984t) and implemented in a slightly modified vari- 
ant h ere). An improvement over previous models of this 
kind (|Kaiserlll983 IValaeeasI [2009h is the good agreement 
at low mass, which is obtained through the constant term 
6o in Eq. ([2"4")l . associated with the normalization of the halo 
bias through Eo. (j2"o]) . Thus, it appears that this constraint 



is sufficient to obtain a good description of the halo bias 
over the whole range cr(M) < 10. Indeed, as explained in 
Sect. 12. 31 since the model is reasonably successful at large 
mass, cr < 1, the integral constraint (|25p ensures that the 
value of the bias over the range 1 < a < 10, which contains 
the other half of the total matter content, has the correct 
magnitude, although the detailed shape has no reason to be 
exact a priori. In particular, at low masses, a > 10, which 
represent a small fraction of the total matter density field 
and do not significantly contribute to the integral (|23|) . the 
normalization (|25[) becomes largely irrelevant (since large 
changes of b(M) would not significantly change the overall 
normalization) so that there is no reason a priori to trust 
our model. Thus, it may happen that at very low masses 
the bias goes to zero, for instance as a power law over M. 
However, this is beyond the reach of current numerical sim- 
ulations and it is not a serious practical problem, for the 
same reason that this only concerns a small fraction of the 
matter content and of the halo population. 

As is well known, the dependence on redshift, at fixed 
<r(M), is quite weak over the range < z < 2.5. There ap- 
pears to be a slight growth of the bias at larger redshift, in 
the regime of rare halos (cr(M) < 0.5). This is most clearly 
seen by the comparison with the dashed line associated with 
the fit from lSheth et al.l ((2001), which is independent of z 
and is identical in the three panels. This small growth is also 
reproduced by our model (1241) . as shown by the solid line. 
This confirms the validity of this approach, and more gen- 
erally of such models that follow IKaiserl (jl984tl . However, 
for practical purposes it is probably sufficient to neglect the 
dependence on redshift in this range. 

It is i nteresting to note that the bias obtained from 
Eq.(8) of ISheth et all (|200lD behaves at large masses as 
b ~ a<5 c /(7g, with a parameter a ~ 0.707 obtained from fits 
to the halo mass function (through the peak-background 
split argument), whereas we have noticed in Sect. l2~B1 that 
Eq.flU]) yields b — (SL*/cr q )(a qiq (s)/a Q:0 (x)). In order to 
explain why both predictions are rather close (especially 
at z = in Fig. [IJ despite these different forms, and with- 
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out introducing such a parameter a in Ea. (|22j) . we have also 
plotted the curves obtained by setting s = x (i.e. neglecting 
halo motions) or Sl* — S c , with the usual value S c = 1.686. 
We can see that the change from 5 C to Sl* makes almost 
no difference for the halo bias (for <5* = 200), while taking 
into account halo motions leads to a significant reduction 
at z — for massive halos. Indeed, since rare halos tend to 
move closer, neglecting halo motions underestimates their 
initial distance and overestimates their correlation. This ef- 
fect is greater at lower redshift, associated with more mas- 
sive objects, and leads to a bias that h a ppens to be very 
close to the prediction from lSheth et al.l ([2001D at z = 0. 

Thus, taking into account halo motions has the same 
effe ct (i.e. reduc i ng th e bias at low z) as the parameter 
a of ISheth et al.l ((2001). Within the peak-background split 
approach, the latter is not a new free parameter for the bias, 
since it is set by the halo mass function itself. However, 
within our approach, we do not need to introduce such a 
parameter, neither for the halo mass function (|29D nor for 
the bias (|22p . We can note that the very weak dependence 
on redshift of the halo bias, in the range < z < 2.5, is 
somewhat misleading, since within this model it arises from 
the partial cancellation between two opposite trends (as the 
"s = x" curve decreases at higher z while the full prediction 
grows to get closer to it). 
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Fig. 2. The mean bias b(x), as a function of distance x, at 
redshift z = for halos defined by the nonlinear density 
threshold <5* = 200. We show our results from Eg . ((30)) for 
the five mass bins 10 13 < M < 3 x lO^/i^Mg, 3 x 10 13 < 
M < W^h-Hlv, 10 14 < M < 3 x lO^h^MQ, 3 x 10 14 < 
and 10 15 < M < 3 x 10 15 h^Mm, from 



M < 10 15 /i- 1 M Q , 
bottom to top. 



In Fig. [T] we have considered the bias at the distance 
x = 50/i _1 Mpc, which is the typical distance where the 
large-scale halo bias is measured in nu merical simulations , 
with box sizes of order SOOft^Mpc (jTinker et al.ll2010ft . 
The bias measured in such simulations is almost scale- 
independent over the ran ge 20 < x < 1 1 0/i~ 1 M pc, as seen 
fo r instance in Fig. 10 of Manera et al.l (|2010l ) and Fig. 17 
of lManera fc Gaztanaeal (|2010D . This agrees with the stan- 
dard large-scale behavior b ~ 5l*/<J 2 . In our model, there is 
an additional prefactor a qtq (s)/ao t o{x), due to halo motions 
and to the fact that high wavenumbers do not contribute 
to halo correlations. Therefore, we show in Fig. [2] the de- 



f£b(M)n(M)dM 



b(x) 



at z = for five mass 



(30) 



We can check that we obtain a weak dependence on scale 
over the range 20 < x < lOOft^Mpc (and we obtain similar 
results at higher redshift). At smaller scales, the upturn for 
the bias of the most massive halos (upper dotted curve) 
is due to the fact that the argument of the exponential in 
Ea.([22|) becomes large so that the linear approximation is 
no longer valid. We shall come back to this point in Sect. [5] 
below. At even smaller scales, of the order of the radius of 
the halos, exclusion constraints that we have not taken into 
account come into play and imply a halo correlation equal 
to —1. It is clear that our model only applies to the larger 
scales, shown in Fig. [21 where halos are well separated. 

At larger scales, the bias is difficult to measure (be- 
cause the two-point correlation functions are quite small) 
an d actually beco mes meaningless. Indeed, as pointed out 
in Valageas (2010) for more general initial conditions, and 
in pre vious studies where halos are identified with densit y 
peaks (|Coleslll989t iLumsden et~aHll989t lDesiacauesll2008h . 
since the halo correlation is not exactly proportional to 
the matter correlation (because of the smoothing at scale q 
in a 2 q (s), of halo motions, and of subleading terms) their 
zero-crossings do not exactly coincide. This implies that the 
ratio b 2 = £(M, x)/S > {x) shows wild oscillations and diver- 
gent peaks at the scale where the matter two-point correla- 
tion changes sign (typically at x ~ 120/i _1 Mpc). Therefore, 
at very large scales, x > 100/i _1 Mpc, it is no longer use- 
ful to introduce the ratio b 2 , and it is best to work with 
the halo and matter correlations themselves (note that this 
conclusion is quite general and not specific to the model 
studied here). 



3. Extension to higher density contrasts 

The results obtained in Sect. [5] applied to halos defined by 
a nonlinear density contrast such that 5* Ss 200, so that 
the large-mass tails could be derived from the spherical col- 
lapse dynamics associated with Eq.([T]) (and we focused on 
the cas e 5* = 200 , whic h is of practical interest). As pointed 
out in IValageasl (|2009f ) , for larger density contrasts shell 
crossing plays a key role and even in the rare-event or large- 
mass limi t EcjJfT]) no longer applies. Indeed, as discussed in 
Valageas (2002), because of a strong radial-orbit instabil- 
ity spherical dynamics is no longer a useful guide. More 
precisely, the properties of the density field in such high- 
density regions are no longer governed by initial configura- 
tions that are spherically symmetric (whence with purely 
radial motions) and one must take into account (infinitesi- 
mally) small deviations from spherical symmetry, that are 
amplified in a non-perturbative manner after collapse. This 
is a significant difficulty for precise and robust modeling 
of halos defined by high density thresholds. Therefore, we 
investigate in this paper a very simple model that tries to 
bypass this problem by relating inner high-density shells to 
outer lower-density radii. 

As we have seen in Sect. [2j in order to obtain a reason- 
ably successful model for the distribution (i.e., the bias) of 
dark matter halos, defined by <5* = 200, it is sufficient to 
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Fig. 3. The halo mass functions at redshift z = 0, for halos defined by a nonlinear densi t y con trast (5* = 200,400,800, 
and 1600. The solid line is our model (1551) . the dashed line is the fit from Sheth fc Tormenl ( 19991 ) . which does not change 
with 5* , and the points are the results from numerical simulations in iTinker et al. ( 20081) . The value of the parameter a 
in Ea. (j3"Tj) is set to a = 2.2. Here M = Ms, is the halo mass enclosed within the radius defined by the density threshold 
<5». 



derive the la rge- mass tail, following the standard ideas of 
iKaiserl (|l984t ) where rare nonlinear objects are identified in 
the initial (linear) density field. Then, normalization con- 
ditions constrain the bias of typical halos (a(M ) i; 1) and 
automatically provide reasonable estimates in this range. 
The advantage of this procedure is to bypass a detailed 
treatment of typical or low-mass halos, which would re- 
qui re taking into account tidal effects and mergings (but 
see iBlanchard et~aTl (|1992[ )). 

In this spirit, in order to extend our model to higher 
density thresholds <5*, we only need to explicitly consider 
rare massive halos. Then, since these halos are isolated and 
relaxed objects, we assume that they can be described by 
a mean mass profile M(< r), with a mean density within 
radius r that behaves as 



P(< 



3M(< r 
4 7r p r 3 



with < a < 3, 



(31) 



where M (< r) is the mass enclosed within radius r. Here we 
could have used a m ore detailed density p rofile, such as the 
usual NFW profile dNavarro et al.lll997|) . as described in 
iHu fc Kravtso v (2003). However, this would also introduce 
a further dependence on mass and redshift through the con- 
centration parameter c(M, z) that parametrizes such pro- 
files. Therefore, in view of the approximations involved in 



our approach we think that such refinements are beyond 
the scope of this study and we prefer to use a single power- 
law index a. This should be sufficient as long as we do 
not consider too wide a range of contrasts (5», over which 
the slope of the density profile shows significant changes. 
Then, from Eq. (l3"l"T) we can relate the mass Mj t , enclosed 
with the radius defined by the nonlinear density contrast 
<5*, to M 2 oo, associated with <5* = 200, by 

/ 201 \ (3 ~ a)/a 
M Sm =M 200 ( — -J . (32) 



With the same one-to-one identification, we can define a 
reduced variable vg, as 

vs m (M s , ) = f 200 (M200 ) , (33) 

where ^200(^200) is given by the second Eq. (|2"5]) and M 2 oo 
is obtained as a function of M$ t through Eq. (|3"2"|) . and the 
large-mass tail of the halo mass function is still given by 
ns t (M$ t ) ~ e~ Vs *l 2 = e _i/ 2oo/ 2 . This one-to-one identifi- 
cation of rare massive halos also implies that their spatial 
distribution remains the same, independently of the choice 
of <5* , so that the large-mass two-point correlation and bias 
are still given by Eqs.([l9]) and (j22"1) . 

bl e ^(M s „x) = bl c . 200 (M 200 ,x). (34) 
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Fig. 4. The halo bias b(M,x) as a function of <j(M), at distance x = 50/i _1 Mpc and redshift z — 0. We show our results 
for halos defined by a density threshold <5* = 400 (left panel), 5* = 800 (middle panel), and <5* = 1600 (right panel). The 
solid l ine is our model (f3"6"|). with a — 2.2 as for the mass functions of Fig. [3] The dashed line is the fit from lSheth et al.l 
(|200l|) . for hal os with 6* = 200, w hich is the same in all panels and in Fig. [T] The points are the results from numerical 
simulations of iTinker et~aT1 (|2010f ). 



Next, in order to describe typical halos we again take 
advantage of normalization constraints. First, to ensure 
that the mass function remains normalized to unity as in 
Eq. (|27[) . we choose the approximation 

n s , (M s ,)dM s , = -f - f{v) — , (35) 

where the reduced function f(v) is still given by Eq. (j2"9")l and 
the reduced variable v is given by Eq. (j35)) . Here the mass 
A/200 is still defined by Eo. (j3"2]) but since there is no longer 
a one-to-one identification (except at very large mass, as 
explained below) A^oo must be seen as an "effective mass" 
rather than the mass within <5* = 200 of the same individual 
object. This prescription automatically provides both the 
right large-mass cutoff, as explained below Eg . (j33j) (pro- 
vided the mean profile pip is correct) and the normaliza- 
tion (|2"T)) . Note that this is a normalization in terms of mass, 
which ensures that by counting such halos we recover the 
mean matter density of the Universe. Thus, a fixed range 
V2] is associated with the same fraction of matter, inde- 
pendently of the choice of <5* , but it is divided into a larger 
number of smaller objects for larger <5*, because of the pref- 
actor p /Ms, in Eq. (135"1) . Indeed, a given v is associated to 
a fixed M200, through the second Eg . (|2"81 . and to a mass 
M s , oc (1 + (5*)-( 3 - Q )/ Q through Eq.([3J). Thus, through 
the normalization constraint and the form (|35[) we include 
some dependence of the halo multiplicity on the choice of 5* 
(as expected, a higher threshold <5* leads to a larger num- 
ber of objects), and the one-to-one identification used in 
Eqs. (|3"2"|) - ([3"4"|) only applies in an asymptotic sense at high 
mass. 

Second, to ensure that the halo bias also remains 
normalized to unity, as in Eq. (j2"5)) . we also add to the 
large-mass term b T . e . ; s, (Ms, , x) a mass-independent term 
bo-,S,(x), as in Eq. (|2"4"l) . Then, thanks to the choice ([351 we 
find that the normalization (|25[) yields bQ-s,(x) — bo-2oo(x). 
Combining with Eo. (l34l) this gives for the total halo bias, 

bs,(M s „x) = 6 2 oo(M 2 oo,a;). (36) 



This automatically satisfies both the expected large-mass 
behavior (where the halos labeled by <5* and 200 are the 
same objects) and the normalization constraint. Since at 
moderate and low mass we no longer have a one-to-one 
identification between halos, as we have seen from the mass 
function (1331) . the equality ([361 should not be seen as the 
result of a strict one-to-one identification, even though it 
happens to give a similar expression. 

We first show in Fig.[3Jthe halo mass functions we obtain 
with our simple model at redshift z = 0. In the following, 
for simplicity we do not write the subscript 5„, so that in 
each panel of Fig. [3] the mass M is actually the mass Ms, 
enclosed within the radius defined by the density threshold 
<5*, as labeled in each plot. 

For <5* > 200 we set the parameter a of Eq. (j3"Tj) to a = 
2.2. Within our approach, this is the only new parameter 
that is needed to obtain the mass functions for arbitrary 
5* from the one at <5* = 200 given by Ea. (j29"]) . We choose a 
value such that the large-mass tail is correctly reproduced 
for 200 < <5* < 1600. This gives a reasonable value for the 
slope of the density profile around the virial radius (i.e. 
p(< r) ~ r~ 2 - 2 ). For larger values of the density threshold, 
which correspond to smaller radii within massive halos, a 
smaller value of a would probably be required in order to 
follow the flattening of the halo profile, down to a ~ — 1. 
Since for practical purposes one usually considers density 
thresholds in the range 170 < <5* < 500 to define dark 
matter halos, we do not investigate further this point. As 
seen in Fig. [31 our simple model is already sufficient to 
reproduce the shift of the large-mass tail over the range 
200 < <5* < 1600. This is most clearly seen by comparison 
with the fixed dashed l ine th at applies to <5* = 200, taken 
from lSheth fc Tormenl (|l999ft . 

For high values of the threshold 8*, especially in the 
lower right panel at 5* = 1600, it appears that our model 
overestimates the mass function for intermediate-mass ha- 
los, a(M ) ~ 1.5. On the other hand, it seems that in numer- 
ical simulations some fraction of matter is lost as we con- 
sider higher thresholds (5* , so that the normalization (|2"T1) no 
longer holds (the integral would be smaller than unity for 
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large 8*). This may have a physical meaning, for instance 
if some regions of space are smooth and show a maximum 
density contrast, in which case they are no longer included 
as we select a higher threshold. A second factor is that 
at some level the distribution of the matter content of the 
Universe over a set of halos is not a very well defined proce- 
dure. More precisely, even though one may always build a 
well-defined algorithm to redistribute particles in a collec- 
tion of halos within numerical simulations, to some degree 
it involves some arbitrariness from a physical point of view. 
This is clear from the fact that one cannot build a parti- 
tion of a 3D box with spheres, and this means that the 
description of the matter distribution in terms of spherical 
halos can only be approximate. This may also lead to a vi- 
olation of the normalization (|27p . A detailed investigation 
of such effects is beyond the scope of this paper and we 
prefer to stick to our simple model and to fixed normaliza- 
tions such as Eq. ([2"Tl) . For practical purposes, the accuracy 
of the model (|35"j) should be sufficient, and appears to be 
quite satisfactory in view of its simplicity. 

We compare in Fig. @] the halo bias predicted by 
our model, Eg .(1361), with results from numerical simula- 
tions (|Tinker et al.l l2010h . for the density thresholds <5* = 
400,800, and 1600. Although the theoretical curve corre- 
sponds to z = 0, in order to increase the statistics we plot 
the numerical data obtained for < z < 2.5, since the de- 
pendence on redshift is quite weak (see Fig. [1]). Of course, 
for the parameter a we use the same value as the one used 
for the mass functions shown in Fig. [3J a = 2.2. 

We can see that as the density threshold <5* increases the 
bias at fixed a(M) grows, especially at large masses. This 
is most clearly seen through the comparison w ith the fixed 
dashe d line, which corresponds to the fit from lSheth et al.l 
(2001) for halos defined by 8* = 200. This moderate depen- 
dence on 8* is also reproduced by our simple model (I361 . 
which shows a reasonable agreement with N-body simula- 
tions over this range, 200 < <5* < 1600. In fact, for a con- 
stant parameter a, as in Fig.0] we can see from Eq. (|36")) that 
as 8* grows the curve b(M) is simply shifted by a uniform 
translation towards the left in the (In M, b) plane, since at 
fixed bias b, whence at fixed "effective mass" M200, we have 
InM = InMaoo + ^ ln[201/(l + 8*)}. Since a(M) is not 
exactly a power law, this horizontal translation is not ex- 
actly uniform in the (ln(l/cr), b) plane of Fig.SJ As noticed 
above, in order to cover a larger range of density thresholds 
(5* , or to improve the accuracy, it would be necessary to use 
a slope a in Eq. (j3Tj) that runs with M and <5*, in order to 
take into account the dependence on mass of the halo pro- 
file and the flattening of the slope at inner radii. However, 
we can see in Fig. 2] that the simple model with a constant 
value, a = 2.2, already provides a good match to numerical 
simulations. 

An advantage of this simple approximation is that it is 
straightforward to satisfy the normalization (|23|) . as shown 
by the simple consequence (|36p . This ensures that our un- 
derlying halo model is self-consistent, in the sense that if 
we desc ribe the matter distribu tion through a standard halo 
model (jCoorav fc ShetbJ I2002T ). but with arbitrary thresh- 
old (5* to define the halos, by integrating over these ha- 
los we recover both the total matter density (through the 
halo mass function) and the matter two-point correlation 
(through the halo bias). This means that, despite the ap- 
proximate nature of such halo descriptions noticed above in 



the discussion of Fig. [3J our model automatically satisfies 
these two integral constraints. 

We think such properties should be taken into account 
in the building of models for the matter distribution. As 
seen in this work, they can serve as a useful guide, which 
can be sufficient to provide reasonable quantitative esti- 
mates over tightly constrained domains. A second benefit 
is that they avoid introducing small inconsistencies, which 
may lead to spurious quantitative discrepancies for quanti- 
ties that would be computed in later steps by integrating 
over the halo populations. 



4. Impact of halo motions 
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Fig. 5. The ratio b s=x /b, as a function of er(M), at redshifts 
2 = 0, 1.25, and 2.5. The bias b is the one obtained in Sect. [5] 
and shown in Fig. [1] whereas b s=x is obtained by making 
the approximation s = x in Eq. 



We now take advantage of the simple analytical model 
presented in Sect. [2] to estimate the impact of halo motions 
on their observed two-point correlation. Thus, focusing on 
halos defined by the nonlinear density threshold 8* = 200, 
we plot in Fig. [5] the ratio b s=x /b, where b is the bias ob- 
tained in Sect. [Hand shown by the solid line in Fig.Q] while 
b s=x is the bias obtained with the approximation s = x (i.e. 
we substitute s — > x in Eq. (l22)) h which was plotted as the 
red dashed line in Fig. [T] Therefore, this ratio measures the 
effect on the bias o f halo motio ns, which are neglected in 
the usual appr oach (lKaiserll984f) o r in the peak-backgr ound 
split method ([Cole fc Kaisej|1989HMo fc Whitdll996l) . 

For rare massive halos, where the asymptotic expres- 
sion (|22[) applies, we have s > x within our complete ap- 
proach, as can be checked from Eo.([2T))). Thus, as could be 
expected, rare massive halos tend to move closer because of 
their mutual gravitational attraction. Although this might 
seem contradictory with large-scale homogeneity (one may 
ask how could all halos move closer to each other ?), this is 
not the case. For instance, let us consider an infinite regular 
3D grid of cell size L, and let us put two halos separated 
by a small distance, I <C L, around each vertex (so that 
the vertex is the center of the pair). Then, if at a later 
time we independently shrink each pair around its vertex, 
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t — > £' with £' < I, we can see that on the mean halos effect of the average halo motion now increases at higher 



have moved closer. Indeed, the distances between different 
pairs have not changed while they have been reduced within 
each pair. However, the system has clearly remained homo- 
geneous on large scales (there is not coherent shrinking of 
all halos positions towards a central point, as one might 
have been afraid of). This simple picture shows how small 
scale motions can have a nonzero effect on the mean, while 
preserving large-scale statistical homogeneity. 



redshift, for a fixed mass M. 

At z = 0, neglecting this effect leads to an overestima- 
tionofthebiasfc(M) by a factor 1.5 at M ~ 3xl0 15 /i _1 M Q , 
that is, cr(M) <~ 0.34. Therefore, it cannot be neglected for 
massive halos. 



5. Impact of exponential nonlinearity 
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Fig. 6. The ratio b s=x /b, as in Fig. [5j but as a function of 
M, at redshifts 2 = 0, 1.25, and 2.5. 



Then, it is clear from Eg. (|2"21 that the predicted bias, or 
the predicted two-point correlation function, of rare mas- 
sive halos is smaller when we take into account the fact 
that s > x than when we neglect this distinction. Indeed, 
cross-correlations, such as (s), of the linear density fluc- 
tuations at scale q over a large distance s are decreasing 
functions of s. This implies that by making the approxima- 
tion "s = x" we overestimate the initial cross-correlation 
between massive halos observed at distance x, and this ex- 
plains why the ratio b s=x /b is larger than unity at large 
masses in Fig. [5] This effect increases at larger masses, 
which have a large bias and show a stronger sensitivity to 
the large-scale correlations. Of course, this agrees with the 
behavior found in Fig. [TJ At small masses the ratio b s=x /b 
becomes slightly smaller than unity. This is a direct conse- 
quence of the normalization constraint (|25[) . which implies 
that the greater value of the bias at high mass, within the 
approximation "s — x" , must be compensated by a smaller 
value at small mass. The effect is rather small in this range 
since massive halos are rare so that their bias does not con- 
tribute much to the normalization (1251). 



We can see in Fig. [5] that the effect of halo motions is 
smaller at higher redshift for a fixed value of a(M), which 
actually corresponds to a smaller mass at higher z. Thus, if 
one only considers typical halos, or a population defined by 
a fixed upper bound on l/a(M) (this roughly corresponds 
to a fixed comoving number density), the average halo mo- 
tion can be neglected at high z. We show in Fig. [6] the same 
ratio b s=x /b at redshifts z = 0, 1.25, and 2.5, but as a func- 
tion of M instead of X/a(M). Then, we can see that the 



Fig. 7. The ratio &£,/&, as a function of cr(M), at redshifts 
2 = 0, 1.25, and 2.5. The bias b is the one obtained in Sect. [2] 
and shown in Fig. [TJ whereas 6l is obtained by using the 
"linear" approximation (|37p in Eq. 
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Fig. 8. The ratio as in Fig. [71 but as a function of M, 
at redshifts z = 0, 1.25, and 2.5. 

At large distance cross-correlations such as cr^ be- 
come very small and one can expand the exponential in 
Eq.(f2"2"j). to obtain the linearized form 



bl c AM,x) 



^k 2 + o-L( s )) °ifi{ x ) 



(37) 
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At zeroth-order ove r <J qq (s ), where s — x, this gives the 
well-known result of iKaiserl (|1984| L b r , e ,.^(M, x) ~ Sl*/^ 2 , 
which shows how the bias of massive halos grows as oc 1 /a 2 , 
in agreement with Fig. [JJ Linearized expressions such as 
(j3"7[) . that is, where one looks for an expression of the halo 
two-point correlation function, or of the halo power spec- 
trum, that is linear over the matter two-point correlation 
or power spectru m, are widely used. In terms of the "lo- 
cal bias model" (iFrv fc Gaztanaeal 119931: iMo et al.l fl997t 
iManera fc Gaztanaga 20101 ). this also corresponds to the 
approximation of "linear bias" , where the halo density field, 
riM (x) , is written as a linear function of the matter density 
field, such as «m(x) = HmW + bi(M)5M(*-)] (where 5m (x) 
is the matter density contrast smoothed on scale M), with 
a constant linear b i as fac tor b±. However, as pointed out 
bv lPolitzer fc Wisd (Il98l . even when the cross-correlation 
a 2 q (s) is small the argument of the exponential can be 
large, because of the factor S^/ a q . For instance, at fixed 
scales q and s the argument grows as 1/c 2 as the amplitude 
of the linear matter power spectrum decreases. Therefore, 
we investigate in this section the impact of the exponential 
non-linearity (|22l) on the halo bias. 

Thus, we show in Fig. [7] the ratio where b is again 
the bias obtained in Sect. [2] and Fig. [TJ while &l is the 
bias obtained with the approximation (1371) . which we sub- 
stitute into Ea. (|!M|) . We still take into account halo mo- 
tions through Ea. (|20p . as well as normalizations to unity. 
Therefore, this ratio measures the effect of the exponential 
non-linearity (|22l) . In agreement with the previous discus- 
sion, we can see in Fig. [7] that by linearizing the expression 
(|2"2"|) we underestimate the halo bias at large masses. Indeed, 
the argument in Ea. ([2"2"j) scales as ~ ^\* (T q,q,s/ (J q K l/°g> 
which grows at large mass. However, this effect is only sig- 
nificant at very large masses and typical objects are not 
affected by this linearization. As shown by Figs. [7J and [FJ 
this effect decreases at higher redshift for a fixed value of 
cr(M), but increases at higher redshift for a fixed mass M . 
Therefore, although for typical halos such a "linearization" 
of the bias is a very good approximation, for very massive 
objects it can lead to a significant underestimation, espe- 
cially at high redshift. As noticed in Fig. [5] (upper dotted 
curve), this effect also becomes more important at smaller 
distance x, as cr\ qB is larger. 

6. Conclusion 

We have described in this paper a simple analytical model 
for the bias of dark matter halos. It extends previous 
works, which focused on halos defined by their virial den- 
sity contrast (6» ~ 200), to the more general case of ha- 
los defined by an arbitrary density threshold in the range 
200 < (5* < 1600. This is coupled to a model for the 
halo mass function that is also generalized to these density 
thresholds, and we have checked that these simple mod- 
els yield a good agreement with numerical simulations at 
redshifts < z < 2.5. 

Our approach also improves over some previous works 
by including the effect of halo motions on their two-point 
correlation function. This arises from the fact that massive 
halos move closer on the mean, because of their mutual 
gravitational attraction, which implies that their mean dis- 
tance at a given time is smaller than the Lagrangian sepa- 
ration of the two regions they originate from in the primor- 



dial density field. We estimate this effect from the linear 
displacement of halos, which provides a simple approxima- 
tion that could be easily used in more complex cases, such 
as non-Gaussian initial conditions. 

Another feature of our model is that it contains no free 
parameters for the bias. More precisely, the only parame- 
ters that appear are related to the halo mass functions. As 
in most other approaches, a few parameters are needed to 
describe the shape of the mass function (of halos defined 
by 5* = 200) at low and intermediate masses, but the large 
mass tail is governed by the cutoff e~ v I 2 without further 
tuning. Then, the extension to halo populations defined by 
larger density contrasts only involves a single parameter, 
a ~ 2.2, which is related to the slope of halo density pro- 
files around the virial radius. Then, no further parameters 
are required to compute the bias of these halo. 

As stressed in this paper, this is possible thanks to the 
use of integral constraints for both the halo mass function 
and bias. This allows us to focus on rare and massive ob- 
jects, where reliable predictions can be derived. Then, the 
behavior of intermediate and small halos, which is beyond 
the reach of analytical approaches because of strong tidal 
effects and mergings, is strongly constrained by these inte- 
gral conditions, and we have shown that this is sufficient to 
build simple and efficient models. Of course, this simplic- 
ity has a cost: we cannot expect a priori high accuracy at 
very low mass, and there is no systematic procedure (such 
as expansions over some parameters) to improve the model 
up to arbitrarily high accuracy. However, we have seen that 
we already obtain a good agreement with numerical simu- 
lations, and this approach should be quite robust. Another 
advantage is that by explicitly taking into account such in- 
tegral constraints we make sure our model is self-consistent, 
and we avoid any risk to introduce spurious discrepancies 
that may arise in integral quantities because of small in- 
consistencies. 

Finally, we have evaluated the quantitative impact on 
halo bias of two common approximations, i) neglecting halo 
motions, and ii) linearizing the halo two-point correlation 
over the matter power spectrum. This could be useful to 
check the range of validity of these approximations. 

It would be interesting to generalize this work to more 
complex cases, such as non-Gaussian initial conditions. 
However, we leave this to future studies. 

Acknowledgements. The author thanks J. L. Tinker for sending his 
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